Isotropisation of flat homogeneous universes 
with scalar fields 



Stephane Fay* 
Jean-Pierre Luminet^ 
Laboratoire Univers et Theories (LUTH), CNRS-UMR 8102 
Observatoire de Paris, F-92195 Meudon Cedex 
France 

7th February 2008 



Abstract 

Starting from an anisotropic flat cosmological model(Bianchi type /), 
we show that conditions leading to isotropisation fall into 3 classes, re- 
spectively 1, 2 ,3. We look for necessary conditions such that a Bianchi 
type I model reaches a stable isotropic state due to the presence of several 
massive scalar fields minimally coupled to the metric with a perfect fluid 
for class 1 isotropisation. The conditions are written in terms of some 
functions I of the scalar flelds. Two types of theories are studied. The 
flrst one deals with scalar tensor theories resulting from extra-dimensions 
compactiflcation, where the Brans-Dicke coupling functions only depend 
on their associated scalar flelds. The second one is related to the presence 
of complex scalar flelds. We give the metric and potential asymptotical 
behaviours originating from class 1 isotropisation. The results depend on 
the domination of the scalar fleld potential compared to the perfect fluid 
energy density. We give explicit examples showing that some hybrid infla- 
tion theories do not lead to isotropy contrary to some high-order theories, 
whereas the most common forms of complex scalar flelds undergo a class 
3 isotropisation, characterised by strong oscillations of the (. functions. 



1 Introduction 



If General Relativity is the best available theory describing our local Universe, it 
has recently become clear P] that the modest amount of matter in the Universe 
(30% of the total energy density) is complemented by a large amount of exotic 
energy (70%). This exotic energy impHes that the Universe is approximately 
spatially fiat, and that its expansion is accelerating. To account for such a 
dynamics, several proposals exist which extend General Relativity to include 
higher order theories j21 El, dissipative fluids^ |^ or massive scalar fields. We 
are going to consider this last type of energy content. 
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Although most of papers only take into account one single scalar field, there 
are many reasons to consider the presence of several ones. Indeed, particle 
physics predicts high-order theories of gravity with extra-dimensions, which can 
be cast into an Einstein form in a 4-spacetime with several scalar fields by 
help of conformal transformations^ EI- In supersymmetry, the adjunction 
of several scalar fields achieves equality between bosonic and fermionic degrees. 
Other reasons may be related to various infiationary mechanisms such as hybrid 
infiation, which needs two scalar fieldsJOjUni: a first one, ip, decreases to a local 
minimum corresponding to a false vacuum. Then the vacuum energy dominates 
and early time infiation begins. During this time, a second scalar field (f) varies 
and when it reaches a threshold value (j)c, a fast variation of if) arises. The two 
fields fit toward some values corresponding to a true vacuum and the end of 
infiation. A last reason could be the presence of complex scalar fields. A scalar 
tensor theory with one complex scalar field ( can be cast into another one with 
two real scalar fields, ip and 0, by help of the transformation ^ = vS™^^™"^' 

From a geometrical point of view, the standard cosmological model lies on 
the assumption that the Universe is perfectly isotropic, homogeneous and thus 
described by the FLRW metrics. However, they are very particular ones among 
the set of all possible metrics and we have to understand why our Universe may 
be described by them. One answer is to assume that it was not so symmetric 
at the beginning of time and that it quickly evolved to an isotropic and ho- 
mogeneous state, as indicated by CMB observations. Moreover the singularity 
approach in FLRW models is far from being generic. Hence, it seems more natu- 
ral to consider that the Universe was born with a more general geometry and has 
evolved toward a FLRW one. One possibility is to leave the isotropy hypothesis, 
keeping only homogeneity. Anisotropic models are described by the nine Bianchi 
models and allow studying how the Universe may tend to isotropy. Their be- 
haviour near the singularity could be shared by inhomogeneous models [TTl [12] 
and one of them admits the fiat FLRW model solution consistent with recent 
CMB observations!^: the fiat Bianchi type / model. 

The goal of this paper is to look for necessary conditions allowing for Bianchi 
type / model isotropisation when two minimally coupled and massive scalar 
fields with a perfect fiuid are considered, and to study the asymptotic dynamics 
of the metric and potential in the neighbourhood of this state. From a technical 
point of view, we will use Hamiltonian ADM formalism giving the field equations 
as a first order differential system. We will rewrite it with normalised variables 
and look for equilibrium stable states corresponding to isotropic ones for the 
Universep4|. Similar techniques have been used in presence of one single mini- 
mally coupled and massive scalar field|l^ and with a perfect fiuid^HI- For both 
cases, isotropic equilibrium points have been found corresponding to power or 
exponential law expansion for the metric functions. Here, we will also examine 
the stability of these results and Wald's cosmological "No Hair" theorem with 
respect to the presence of additional scalar fields. Intuitively, one could think 
that nothing should change and that the generalisation implying several scalar 
fields should be straightforward. However, we will see that it is not always the 
case and depends on the form of the scalar-tensor theory with respect to these 
fields. 

The paper is organized as follows. In section 12 we derive the Hamiltonian 
equations and rewrite them with normalized variables. In sectional we explain 
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the assumptions we will use to study these equations. In section^lwe determine 
the equilibrium points, the monotonic functions and the asymptotic behaviour 
of the metric near equilibrium. We summarize and discuss our results in section 
El In sectional some explicit applications are performed and we conclude in 
section 13 



2 Field equations 

In this section we calculate the field equations. The metric of the Bianchi type 
/ model is: 

ds^ = -{N^ - N,N')dn'^ + 2N,dnuj' + Rlg.juj'uj^ (1) 

where the are the 1-forms defining the homogeneous Bianchi type / model. 
The gij are the metric functions, N and Ni respectively the lapse and shift 
functions. The relation between the proper time t and the time SI is dt^ — (iV^ — 
NiN^)dW^. In what follows we rewrite the metric functions as gij — e^2Si+2ftj 
and use the Misner parameterisation^^ defined as: 

= diag{l3+ + V3/3_ , (3+ - \/3/3_ , -2/3+) (2) 

pI = 2^^l - 2/ii,6ln\ (3) 
6py = diag{p+ + VSp^ , p+ - Vip- , -2p+) (4) 
the pij being the conjugate momenta of /3y . Hence the metric is cast into: 

ds^ = -(iV^ - N^N')dn^ + 2N,dnuj' + i?^e"2n+2ft,^»^i 
The most general form for the action is: 

S = (16^)"i J[R- (3/2 + c^)0^^0,^0-2 _ (3/2 + ^)^-f'^^^^-^ _ U 

+ 16TTc^L^]V^d'^x (6) 



where a prime denotes ordinary derivation, (j) and ^p are two scalar fields and 
their Brans-Dicke coupling functions with the metric are Lu{(f),tp) and fi{(j),ip). 
U{(l), ip) is the potential and the Lagrangian of a perfect fiuid with an equa- 
tion of state p = (7 — l)p- We will consider the interval 7 € [1,2] in which 
7 = 1 stands for a dust fiuid, 7 = 4/3 for a radiative fiuid. Vacuum energy 
corresponds to 7 = and is equivalent to the presence of a cosmological con- 
stant which we will study. Defining the 3-volume V hy V = e~^^ and using the 
energy impulsion conservation law for the perfect fiuid, = 0, we get for its 
energy density p — . 

Hamiltonian ADM formalism^HI ^| needs to rewrite the action under the fol- 
lowing form: 

S = (16^)-i I K^-^ + 7:^^ + ^''^- - N,C^)d^x (7) 



3 



TTij , TT^ and TT^ are respectively the metric functions and scalar fields conjugate 
momenta. In the action JJjl, N and Ni play the role of Lagrange multipli- 
ers and C° and C" are respectively the superhamiltonian and supermomenta. 
Considering Jsj and we deduce that: 





(8) 
(9) 



where | means covariant derivative on a = const} surface. The variation of 
the action with respect to Lagrange multipliers leads to the constraints C" = 
and = 0. Using Misner parameterisation and the above definition of gij, 
we redefine the action as 5* = J p+dp+ + p-dP- + p^dcj) + p^dip — HdQ 
with p^ = TTTT^, p^ = nn^ and H = 2mr^, the ADM Hamiltonian. Then, the 
constraint = 0, yields for H: 
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From HIOII . we derive the Hamiltonian equations: 
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The dot means time derivative with respect to fl. We choose the shift functions 
such that iV* = and we find that the lapse function is related to the metric 
and the Hamiltonian by the relation d^/dO, ~ — l/27r^7V ^^l- Thus, it comes: 

N = (18) 

The relation between 17 and the proper time t is then dt — —NdVl. We wish to 
rewrite the Hamiltonian equations with normalised variables. The Hamiltonian 
lfTn|. which also stands as a constraint equation, leads to the following choice: 

X = H-^ (19) 

y ^ ^fPmjn-^ (20) 
z =p^0(3 + 2t^)-i/2i/-i (21) 

w =p^^V'(3 + 2Ai)"^/2i/~^ (22) 

It implies that /7>0,3 + 2(jj>0 and 3 + 2/i > so that the variables be real 
and the weak energy condition be satisfied. In addition, we define a variable 
depending on the above ones: 

Some of these variables may be physically interpreted, x is proportional to the 
shear parameter S defined in Ql] and /c^ to the density parameter 57™ of the 
perfect fiuid. We introduce them in the constraint ljl()|l and get: 

V^x^ + i?2y2 + I2z2 + 12u|2 + = 1 (23) 

with = P+ + P- and i?^ — 247r^i?Q. This equation shows that the new 
variables {x,y,z,w) are normaHsed. Then, we rewrite the field equations as: 

x = 3R^y^x-3/2{-f~2)k^x (24) 

y = y{6e^,z + 6e^,w + 3R^y^ - 3) - 3/2(7 - 2)k^y (25) 
i = y^R^{3z - 1/2^0,) + 12w{we^^ - zi^^) - 3/2(7 - 2)k^z (26) 
w = y^R^{3w - 1/2%) + 12z(z% - wi^^) - 3/2(7 - 2)k^w (27) 
where we have defined the folowing functions of the scalar fields (j) and ip 

e^, = 0[/^(7-i(3 + 2co)-i/2 
= ^U^U'^i3 + 2fi)-^^^ 
£4., = 0/i^(3 + 2^)-i(3 + 2cj)-i/2 
e^, = ^Puj^{3 + 2uj)-\3 + 2fi)-'^/^ 



Remark that these equations are unchanged under the transformation x ^ —x 
and/or y ^ —y. Hence, we can limit our study to positive x or y. Moreover, 
some first degree equations for the scalar fields (which are not normalized) may 
be written as: 
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Hence, the nine Hamiltonian equations are rewritten with the six equations 
II24I29|I . This reduction in the number of equations comes from the fact that 
equations ifnji imply p± consts and thus (3+ oc (3-. Hence, it stays 9-3=6 
equations to solve. 



3 Assumptions 

In this section we describe the assumptions we will use to study the above equa- 
tions system. They concern the dependence of the Brans-Dicke functions and 
the potential with respect to the scalar fields, the type of equilibrium isotropic 
states we will consider and how fast it is approached by this system. 

We will study the two following classes of theories from the Bianchi type / 
isotropisation viewpoint: 

• For the first one, uj and /i will respectively depend on (j> and only, 
i.e. -^02 = = whereas U will depend on both scalar fields. It 
means that the coupling between them only appears via the potential. As 
pointed in the introduction this type of theories may be obtained when 
one studies the hybrid infiationjEl or as the outcome of extra-dimensions 
compactificationj^. The theories of jH] and [Jj are commented and studied 
from the isotropisation point of view in respectively the sections 16.11 and 

• For the second one, U and will depend on tp only whereas uj will contain 
both scalar fields. Then, we will have = = 0- This type of theories 
is obtained when one casts a Lagrangian with one complex scalar field 
into another one with two real scalar fields. Complex scalar fields have 
been studied in (20] where the scalar fields quantization is considered, in 
|2H to study the formation of topological defects and in for the Bose- 
Einstein condensate. We have analysed each theory of these papers from 
the isotropisation point of view in respectively the sections 16.31 16.41 and 

Looking at the field equations, we have identified three types of isotropic equi- 
librium states (all characterised by a; ^ when 17 ^ — as as we will show it 
below) that we have classified into three isotropisation classes: 

1. Class 1 is such as all the variables but not necessarily the scalar fields reach 
equilibrium with y ^ 0. Mathematically, it is the only one which allows to 
fully determine the asymptotical behaviours of the metric functions and 
potential in the vicinity of the isotropy. 

2. Class 2 is such as all the variables but not necessarily the scalar fields 
reach equihbrium with y = 0. It is generally not possible to determine the 
asymptotical state of the system near isotropy because of y vanishing. If 
it is technically possible, the study of the general properties of this class 
will be the subject of future work. 
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3. Class 3 is such as at least x± reach equilibrium but not necessarily the 
other variables. If one of them behaves in this way, since it has to be 
bounded as 17 ^ — oo, it would mean that it should be oscillating but not 
damped and then its first derivatives should oscillate around 0. It can 
never happen if one of the £ diverges monotonically or with sufficiently 
small oscillations since then, at least two of the derivatives y, z or w 
will keep the same sign and thus will not be oscillating. It does not 
arise if the £ tend to some constants which is confirmed by numerical 
simulations. However, a partial equilibrium may occur for sufficiently 
oscillating £ which then allow an oscillation of the sign of the (y, z, w) first 
derivatives although x± tend to zero. 

In this paper we will only study the first type of isotropic equihbrium state for 
the following reasons. Mathematically it is the only one allowing to determine 
completely the asymptotical behaviours of the metric functions and potential 
in the vicinity of the isotropy. Physically, near the isotropic state, either one of 
the scalar fields energy densities will be neghgible with respect to the other or 
they will both behave in the same way. Let us assume without loss of generality 
that near the isotropic state the dominant scalar field energy density be the one 
of 4>. y is then proportional to (p^ — p^) / Huhhle^ where Hubble is the Hubble 
function. Defining the scalar field parameter as Q.^ oc pcf,/ Huhhle^ , the class 2 is 
thus such as ^ or whereas the class 3 should be such as fi^ does 

not reach equilibrium. The class 1 is thus the only one such as asymptotically 
tends to a non vanishing constant. WMAP observations indeed shows 
that today il^ — 0.73 and Pct,/P4> < —0.78. 

As a last assumption, we will suppose that the Universe approaches suf- 
ficiently fastly its isotropic state. This is a reasonable assumption since 
the Universe was already very isotropic at the CMB time and it will allow us 
to recover classical behaviours for the metric functions in the vicinity of the 
isotropic state, such as power and exponential laws of the proper time. All the 
asymptotical behaviours we will determine will be concerned by this assump- 
tion. Mathematically, it means that on one hand a function / of the scalar 
fields and on the other hand the variables {y,z,w,k) will have to tend suffi- 
ciently fastly to their equilibrium values such as their variations in the vicinity 
of the equilibrium may be neglected. 

The form of the function /((/>, -0) will be related on the presence or not of a 
perfect fiuid and the dependence of the Brans-Dicke coupHng functions and 
potential with respect to the scalar field. If in the neighbourhood of the equi- 
librium, / tends to a constant equilibrium value /o, vanishing or not and such 
a.s f ^ fo + Sf with Sf << /o, we will assume that / fdfl foft -\- fi, fi 
being an integration constant. It will be equal to the constant /i if /o = 0. This 
assumption could be easily raised by keeping the integral but then our results 
will not be on a closed form and not easily physically interpretable. However it 
is mathematically feasible. 

The same kind of assumptions will be made for the variables {y,z,w,k) with 
respect to {5y, Sz, Sw, 5k) but they can not be raised so easily. A perturba- 
tive analysis would be probably necessary and could depend on the particular 
form of the Brans-Dicke functions and potential with respect to the scalar fields 
whereas we wish to keep these functions undetermined. 
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The above assumptions are illustrated by an example in section B.l.ll in the part 
"Asymptotic behaviours". In the section El where we will apply our results to 
some scalar-tensor theories, they will be systematically checked. 

We will examine each of the above defined classes of scalar tensor theories, 
firstly without a perfect fiuid (fc = strictly) and secondly with it (fc ^ or 
A: ^ 0). Equations II28I29|I will serve to establish the scalar fields asymptotic 
behaviours. 

4 Study of the equilibrium states 

We are going to look for the equilibrium points representing an asymptotically 
isotropic Universe for the two classes of scalar tensor theories defined respec- 
tively by = = and i^^ = =0 and such as all the variables (x, y, z, w) 
reach equilibrium with y 0. 

In their famous paper [24|, Collins and Hawking defined the isotropy as ^ 
— oo, in the following way 

• Let Taf3 be the energy-momemtum tensor: T^^ > and 

rpOi 

ytru represents a mean velocity of the matter compared to surfaces of 
homogeneity. If this quantity did not tend to zero, the Universe would 
not appear homogeneous and isotropic. 

• Let be o-y = {de^ /dt)k(i{e^'^)j')k and — aija^^: ^j^^ 0, i.e. the 
shear parameter, proportional to x variables disappears. This condition 
says that the anisotropy measured locally through the Hubble parameter 
Hq tends to zero. 

• P tends to a constant /3o 

This condition is justified by the fact that the anisotropy measured in 
the CMB is to some extent a measurement of the change of the matrix /3 
between time when radiation was emitted and time when it was observed. 
If (3 did not tend to a constant, one would expect large quantities of 
anisotropies in some directions. 

Hence, in our calculations, we will look for the equilibrium states respecting 
the second point, i.e. isotropisation occurs when a: ^ as ^ — oo. It is 
thus a stable state arising for a diverging value of t. These two limits do not 
depend on each other and their consistency will have to be checked. We will 
see that the third point will be always respected since /3± will always disappear 
exponentially. The first point is also respected since we consider a diagonal 
tensor T"'^ and positive energy densities for the perfect fiuid and scalar field. 
The system of first order equations II24I27|I is not totally autonomous since , 
(4,2, iipi and some functions of (f) and ip. To make it fully autonomous, 

we have to consider the two additional first order equations Ij28l29ll for (p and ip. 
Since the scalar fields do not appear in the constraint, they do not need to be 
bounded whatever the isotropisation class. Hence, looking for stable isotropic 
states for class 1 isotropisation only consists in finding the values of (x, y, z, w) 
depending on the scalar fields such as (i, y, z, w) — (0, 0, 0, 0). The equihbrium 
values of z and w will be introduced in equations Ij28|l and Ij29|l to respectively 
get 4> and ^p asymptotic behaviours. 
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4.1 Without a perfect fluid 
4.1.1 = = 

In this subsection and the following ones, we first look for equilibrium points 
corresponding to isotropic stable states, i.e. such as a; = 0. Then we search for 
monotonic functions and finally calculate the asymptotic behaviours of some 
important quantities in the neighbourhood of these points. 

Calculus of the equilibrium points. 
We find two equilibrium points: 

{x,v,z,w) = (0,±(3-^2^ -4)i/2(V3i?)"\£^,/6,£^,/6) (30) 

and a set of points defined hy y ~ Q. The two first ones respect the constraint 
and are real if + tends to a constant smaller than 3. Both and 
have to tend to a constant such as z and w vanish. We will show below that they 
are in agreement with a negatively diverging value of Vl. We do not look after 
the set of points defined by y = since it concerns the class 2 isotropisation 
which we will not study in this work. 

Monotonic functions. 

By examining equation Ij24|l . we deduce that a; is a monotonic function of 0: 
if a; < (x > 0) initially, it will keep the same sign and will be decreasing 
(increasing). Thus, if initially the Hamiltonian is positive, from the definition 
itTsll of the lapse function TV and the relation between the proper time t and SI, 
we derive that — oo corresponds to late times epoch. In the same way, if 

initially z < ^0^/6 (z > f^^/G), z will be monotonically decreasing(increasing). 
The same conclusion arises for w related on the value l-^^j^. 
As in the case of a single scalar fieldJHI, x being monotonic and with a constant 
sign, we are able to show that the metric functions may have one extremum at 
most. Indeed, their derivatives write as dg^jdl = —2N~^e~^^~^'^^^^ — 1) and 
Pij is a linear combination of /3± which depends on the monotonic function x. 
Consequently, there exists only one value for x such as dgij/dt vanishes. In a 
general way, if we consider the two Brans-Dicke coupling functions as depending 
on both scalar fields (j) and -0 (i.e. ii ^ whatever i = (j)2,ipi,'ip2), z and w 
are not necessarily monotonic but it is always the case for x. Thus, whatever 
the dependence ofuj, ^ and U on (j> and ip, the metric functions will always have 
one extremum at most. 

All these elements show that there is no periodic or homoclinic orbit in the 
phase space {x, y, z, w). 

Asymptotic behaviours 

As explained in the section we define the function / = + such as 
when U, — oo, we assume that / fdfl (^^^ + ^^JSl + /i. Then, when 
diverges and after having replaced y by its equilibrium value neglecting 
its variation Sy in the vicinity of the equilibrium, we deduce from lf2ill that 

^ 4>i ^i' ^ xqc^ ''■I 'I'l' . It shows that the equilibrium points 
reality condition is in agreement with the vanishing of x when SI —oo. Intro- 
ducing this asymptotic expression for x in the lapse function Ijl8|l , it is possible 
to calculate the asymptotic form of as a function of the proper time t, i.e. 
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the metric functions attractor. Its local or global nature can not be determined 
unless we specify uj, fj, or U. When tends to a non vanishing constant, 

e ^ V -t-i^ i'l' . When it vanishes, e~ tends to an exponential oft. We also 
evaluate the asymptotic forms of (f) and ip by rewriting the equations Ij28|l and 
(I29II near the equilibrium. We get two differential equations whose asymptotic 
solutions take the same forms as those of and V when ft ^ —oo: 

1 20^^, ^3^^ 



^ (3 + 2uj)U 
■ _ 2i.^U, 

^-{i + 2^,)u ^^^^ 

Since U = U^(f> + U-^ip and applying our assumption on /, we determine that 
near equilibrium: 

u ^exp[2{el^+eljn] (33) 

This shows that if ^^^^^ ^ vanishing constant, U asymptotically 

vanishes as If it vanishes, the potential tends to some constant. Note that 
the special case for which + ^3 impHes y ^ and thus belongs to class 
2 isotropisation which will not be studied in this work. Approach of equilibrium 
is represented by a phase portrait diagram on figure Q 

4.1.2 e^^=e^^^o 

We proceed as in the previous subsection. 
Calculus of the equilibrium points. 

We find the following equilibrium points, Ei and E2, which might correspond 
to some isotropic stable states: 

El = (0,±(1- 4/3)1/2^-1^ o,^v.i/6) 
E2 = i0,±[2i^,ii^,+2i^,)-^]'^^ R-\ 

+ 21^,1^^ - 3)1/2 J2V3(£^, + 2£^J 



They both check the constraint equation. The first one will be real and bounded 
if < 3 and tends to a constant. The second one needs that i^p^i^iii +2^02)"^ 
tends to a positive constant, + 2iTp^) > 3 and + 2£^2 ^ 0. Remark 

that for E2, i^-i and i-^^ may be unbounded. A third set of equilibrium points 
is (y, z) = (0, 0) but we discard it for the same reasons as in the previous sub- 
section. 

Monotonic functions. 

As written in subsection l4.1.1[ a; is a monotonic function of fl and fl{t) a mono- 
tonic function of the proper time whose limit fl —00 corresponds to late 
time epoch when the Hamiltonian is initially positive. If > {^■4>2 < 0) 
and w > £^i/6 {w < £^^/6), w is an increasing (decreasing) function of SI. If 
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Figure 1: Case lA - Equilibrium point approach when no perfect fluid is present 
and {L^-^, L^2, L^-^, L^2, R,p) — (0.23,0,1.58,0,2,1). The point is located at 
{x,y,z,w) = (0,0.19,0.04,0.26). 



moreover > (i?^^ < 0), w is positive (negative) and keeps a constant sign. 
Asymptotic behaviour in the neighbourhood of Ei 

Here we define / = and write that in Q — oo, J i'^^dQ i'^^i} + fi. 
Then, from (12611 we get: 

z f.i^-'^l^)^-'^ J 1^4.1^2'^^ (34) 

Indeed must tend to a constant but l^^ may diverge. It is why an in- 
tegral of appears in this last expression. It shows that we must have 
(3 — t'^^)^ ~ 2 / l^p^t^r^dO. — oo when Vl — oo such as z vanishes. Moreover, 
considering equation (l27jl where a z^l^^ term is present, we deduce that z has 
to vanish sufficiently fast to allow w equilibrium, i.e. z^t^^ 0. When the 
condition for z vanishing is respected, zz = {3 — (■'^-^)z'^ — 2^^j^^2Z^ and 
we deduce that z^t^^ ^ is always true as long as (obviously) does not 
diverge or/and does not vanish. Otherwise, nothing can be deduced from 
zz vanishing. 

The variable x behaves as a;oe^ '''i'^ and vanishes as ^ —00 when real- 
ity condition for the equilibrium points is respected. As previously, using the 
expression for the lapse function and the relation dt = —NdQ, we get as 
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a function of the proper time near isotropy. If tends to a non vanishing 

constant, tends to t "^i . If vanishes, tends to an exponential of the 
proper time. In the same way as in subsection l4.1.1[ we calculate the differential 
equations whose solutions asymptotically correspond to the forms of (f) and %p 
when r2 — s- —oo: 



Since U — U^^, it comes that: 



iJcxe'^^i" (37) 



Thus, the potential behaves as i ^ when tends to a non vanishing constant 
or as a constant when tends to zero. Again if i'^^ ^ 3, y vanishes and thus 
isotropisation is of class 2. We thus exclude this value from our study. 

Asymptotic behaviour in the neighbourhood of E2 

We define / = £^i(€^i + 2£^2)~^ ^-nd write that in fl ^ -co, J ^^^(^^i + 
2£^2)~^dfl {i^-^ + 2^02)"^^ + /i- With this assumption we calculate that 

near E2, x behaves as a;oe^[^^"*2(^^i+2^"*2)"']f^. Since 2i^^{i^^ +2€^J-i tends to 
a positive constant, it thus vanishes as ^ -co. Using the expression Ijl8|l for 
the lapse function, we calculate that when the quantity 1 — 2^^2 (£^j +2£^2)^^ — 
^tpii^tpi + 2i^r^)~^ tends to a non vanishing constant, ^H^^+2i^^){M^^) _ 

From reality condition for the point E2, we deduce that this power of t is 
positive. Hence, this last expression is increasing when the proper time t di- 
verges in accordance with the growth of when —00. If the quantity 
(^i>i{£tl)i + 2^02)"^ vanishes, the metric functions tend to an exponential of the 
proper time. From (I12II and Ijl3|l we derive the differential equations allowing (j) 
and ^ to get asymptotic forms: 

J, ,/-3C/2(3 + 2Ai)(3 + 2uj) + ^P^U^ [[/(3 + 2u;)]^ 

^ = -2\/3^ ^ ^— — — ^ (38) 

[C/(3 + 2u;)]^ ^ > 

■ 6[/(3 + 2cj) , , 

[U['d + 2uj)\^ 

This last equation easily integrates to give ?7(3 + 2(jj) — e^^^^^^^"), fig being 
an integration constant. Since we have U = U^tp, we deduce from (l39)l that 
jj _ g6£^j(£^j+2^^2) ^. Consequently, when + 2£^2)~^ tends to a non 

vanishing constant, the potential vanishes as t~'^ . If it vanishes, the potential 
tends to a non vanishing constant. 

Approaches of both equilibrium points are represented by phase portrait dia- 
grams on figures 12 and 

4.2 With a perfect fluid 

There are two types of equilibrium points when we take into account a perfect 
fiuid depending if k, or equivalently the density parameter of the perfect fiuid. 
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Figure 2: Case 2 A - First equilibrium point approach when no perfect fluid is 
present and (L^^ , L^^ , L^^ , L^^, = (0,0,0.53,1,2,1). The point is located 
at {x,y,z,w) = (0,0.47,0,0.09). 



tends to a non vanishing or vanishing constant. The first type is studied in the 
two next subsections and the second one in the third subsection. 

4.2.1 e^, = = 

Calculus of the equilibrium points. 

In the annexe 2, we look for the zeros of Ij24l27ll and introduce them in the 
constraint to determine k. The only ones in agreement with isotropy are: 

i?4,5 = (O,±l/2V3i?-i[7(2-7)(^|i+4i)-^]'^',l/47V(^0i+4i)"'> 

with fc^ ^ 1 — 2{e'^^1e'^ ) that is real and non vanishing if + £'^^ > 8/27. 

We will show in this subsection that k tends to non vanishing constant and it 
is why we exclude the special value i'^-^ + i'^-^ 8/27. Since we study the class 
1 isotropisation, y is non vanishing and then + can not diverge. Hence, 
from the forms of -£^4^5 points, we deduce that respectively the sum i'^^ + 
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Figure 3: Case 2 A - Second equilibrium point approach when no perfect fluid 
is present and (i^^, L^j, L^^, L^^, — (0, 0, 4, 1, 2, 1). Let us note how this 
approach is different from the first equilibrium point, x and y undergo damped 
oscillations when they approach their equilibrium values. The point is located 
at {x,y,z,w) = (0,0.29,0.22,0.08). 



and its individual components have to tend to some constants. 

Monotonic functions 
Equation Ij24|l shows that x is a monotonic function with a constant sign what- 
ever the values of , i-ip^ ) ^^i ^^'^ ^ipi • Consequently the lapse function also 
has a constant sign and f2 is a monotonically decreasing function of t if initially 
H > 0, tending to —00 for late times. If z > ^^i, z is a monotonically increasing 
function. However, nothing can be deduced when z < because of the perfect 
fluid presence. The same reasoning holds for w with respect to £^i. Hence, it 
seems that no periodic or homoclinic orbit exists. 

Asymptotic behaviours 
Here, there is no need to make any assumptions related to a function / as defined 
in the subsection El Using Ij28l29ll . the scalar fields asymptotic behaviours are 
defined by the asymptotic solutions of the following systems: 

(3 + 2^i)ct>^UU^ 
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ijj — 37- 



(3 + 2uo)(j)i!UU^ 



(41) 



Linearising II24II in the neighbourhood of i?4,5, we find that asymptotically x 
e~2(T-2)0 ^jjj vanishes as ^ —00 for the considered interval of 7. Then, from 

— o — 

the lapse function TV, it comes that e ^ t^-i . We deduce from y definition and 
the fact that this variable does not vanish near equihbrium that V — > cx V~'^ . 
In the same way we deduce from k definition, using the asymptotic expressions 
for Uii) and that it tends to a non vanishing constant. It is thus the same 
for the density parameter of the perfect fiuid. Approach to equilibrium is 
represented by the phase portrait diagram on figure 21 
Q 





Figure 4: Case IB - Equilibrium point approach when a perfect fiuid with 
7 = 1 is present and (L^^^L^^^L^^^L^^^R^'p^k) — (0.23,0,1.58,0,2,1,0.41). 
The point is located at (a;,?/,z,w) = (0,0.27,0.022,0.15). 



4.2.2 £^^=£^^=0 

Calculus of the equihbrium points. 

We proceed as in the previous section. The details for the equilibrium points 
calculus are given in the annexe 2. We find that the only ones corresponding to 
an isotropic state are: 



i?2,3 = (0,±l/2i?-i^TV37(2-7),0, 1/47C) (42) 
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with fc^ ^ 1 — 3/27^?^^^ and thus require i'^-^^ > 3/2j. We will show below that k 
is non vanishing and it is why we exclude the value 3/27. Moreover, since 

we consider a class 1 isotropisation, y can not tend to zero and is bounded. 
Hence, equilibrium is reached only if £^1 tends to a constant such as y and i 
vanish. 



Monotonic functions 
As already noted in the previous section, x is a monotonic function of fl and has 
a constant sign. Consequently, is a monotonic decreasing function of t if ini- 
tially the Hamiltonian is positive and fl ^ —00 corresponds to late times epochs. 

Asymptotic behaviours 
Once again, there is no need to make any assumptions related to a function / as 
defined in the subsectionEl Surprisingly, equilibrium points £^2,3 have the same 
form as in the presence of a single scalar field -0 in ESI- If we consider a La- 
grangian with a single complex scalar field C and cast it into another Lagrangian 
with 2 real scalar fields and ip, £'2,3 would only depend on its amplitude ip 
and not on its phase </>. 

Again, we find that asymptotically x e"*^'''^^^^ and thus e^^^ t~ in- 
dependently on the scalar fields behaviours. Hence considering the definition 
of y, the potential vanishes as cc V~'^. As in the previous section, one can 
show using these asymptotical behaviours for the metric functions and potential 
that k tends to a non vanishing constant. The differential equation giving "0 in 
f2 —^ —00 may be written: 

^ = 377^ (43) 

To determine a similar equation for cj), we need to know z asymptotic behaviour. 
We find z e^l-*"^ 7/2)0 7 Jf^2f^,i<io] ^^^^^ vanishing when — > —00 

if (1 — 7/2)r2 — 7/ i,p2i^ldn —> —00. Then, the differential equation giving the 
asymptotic form for is: 



120 3[(l-7/2)n-7 / 



^o|_^i-7/z;si-7 J t^2t^irfOj ^^^^ 



V3-K2W 



00 being an integration constant. As in the absence of a perfect fiuid, z has 
to vanish sufficiently fast such that w reaches equilibrium, i.e. we must have 
z^i^2 0. This condition is always satisfied as long as the one allowing the 



vanishing of z is respected, since we have then zz = 3{l — j/2)z^ — 3^£^2£^iZ^ 
and does not diverge. Approach to equilibrium is represented by the phase 
portrait diagram on figure 



z 



4.2.3 The case fc ^ 

As shown above, the limit fc ^ disagrees with the isotropic equilibrium states 
defined for the equilibrium points of subsections 14.2.114.2"^ However, it is al- 
ways possible to assume fc ^ in the field equations and then to solve them. 
We thus recover the equilibrium points obtained in the absence of a perfect 
fiuid. The asymptotic behaviours of x and of the metric functions are the same 
as in section ITU However, the conditions for isotropisation are modified since 
now k has to vanish asymptotically, thus representing an additional constraint. 
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Figure 5: Case 2B - Equilibrium point approach when a perfect fluid is present 
anci {L^-^, L^2, L^-^, L^2, R,p,k) = (0,0,1.53,0.23,2,1,0.60). The point is lo- 
cated at {x, y, z, w) = (0, 0.28, 0, 0.16). 



To flnd it, we rewrite k as fe^e^^^ Then, in the case iff,^ = = 

where x e"^"^*!"^*! , k will vanish only if i'^^ + l"^^ < 3/27 which is consis- 
tent, although more restrictive, with reality condition of the equihbrium points 
when no perfect fluid is present since 7 G [1,2]. Hence, + > 3/27 is 
a necessary condition for isotropisation to occur with fc ^ toward the equi- 
librium points of subsection 14.2. 11 whereas + < 3/27 is a necessary 
condition for isotropisation to occur with k ^ toward the equihbrium points 
of subsection 14.1.11 The same reasoning may be followed concerning the case 
^(fi — ^<t>2 = 0- For the El point, k vanishes only if < 3/27 and for the E2 
point, if 2^0,(f0,+2£,J-i > 1 - 7/2 with 1 - 7/2 e [0, 1/2]. 
Since k — Sy^U^^V^'^ and we consider the class 1 isotropisation such as y ^ 0, 
k vanishing implies U >> and thus, in the Lagrangian fleld equations, the 
potential will dominate the perfect fluid energy density term. 

The results of this last subsection are more accurate and extended than 
those we had found in ^gj. In this last paper, we had considered different cases 
depending on the behaviour of U with respect to and then deduced this 
of k. In the present paper, the opposite reasoning is made and seems to give 
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better results. In particular in 121, we had not detected that conditions for 
isotropisation were changed when fc ^ with respect to the no perfect fluid 



4.3 Technical results summary 

In this subsection, we summarise our technical results. We got the hamiltonian 
field equations for a Bianchi type I Universe filled with a perfect fiuid and two 
scalar fields defined by £ci>2 — ^ip2 — ^^'^ ^<l>i — ^>p2 — 0- We rewrote them with 
normalised variables and looked for the stable isotropic states defined such as 
the shear disappears when the Universe expands, i.e x ^ when SI — oo. We 
then found the equihbrium points summarising in tableland depending on the 
asymptotic behaviour of k or equivalently the perfect fiuid density parameter 



^02 — ^ll>1 — ^ 



, _^ [2U2 (<;Y.i+2^fe)'']''" 

= (U, ± , 



or 

f^,n 



^ ' VSR ' 6 ' 6 



n,n — ^ const 
with 
const ^ 



1/2 



(0,±275J;7v/37(2-7),0,4j^) 



Table 1: The {x,y,z,w) equihbrium points representing stable isotropic states 
for the Universe 



We then found some asymptotical necessary conditions for isotropy depending 
on inequality and limits written with respect to the functions £ of the scalar 
fields. Prom the viewpoint of asymptotical behaviours, it comes asymptotically 
that either the potential vanishes as or tends to a constant. In this last 
case, the Universe tends to a De Sitter one whereas when tends to a non 

2 

vanishing constant, the metric functions behave as t"-' as in the absence of any 
scalar fields. In the other cases, they behave as some powers of the proper time, 
the powers beeing some constants defined as asymptotical limits of some scalar 
fields functions summarised in table El 



£(p2 — ^l/)2 — 




^2 _|_ f2 


E2 : {e^, + 2e^„){3e^j-^ 



Table 2: When the potential desappears and flm does not tend to a constant, 
the metric functions behave as some power laws of the proper time in the neigh- 
bourhood of the isotropic state. These powers are summarised in this table for 
the corresponding isotropic equilibrium point. 
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5 Discussion 



After the tedious computation of previous sections, we now summarize and dis- 
cuss our results. We have classified isotropisation process into three classes and 
looked for necessary conditions allowing for class 1 isotropisation of Bianchi type 
/ model when two minimally and massive scalar fields with a perfect fiuid are 
considered. The class 1 is such as all the variables {x,y,z,w) reach equilibrium 
with y 0. We have assumed that the potential is positive, the scalar field re- 
spects the weak energy condition and the Universe isotropises sufficiently fastly 
such as we could neglect the variation of /((/), "0) and {y, z, w, k) in the vicinity 
of the equilibrium. 

The first necessary conditions we have found for isotropisation stem from the 
definition of isotropy: it will only happen for a vanishing shear (x 0) and an 
eternally expanding Universe {Vt — oo), thus implying that an isotropic state 
is always stable. Moreover, it will correspond to late time isotropisation if the 
Hamiltonian is initially positive. Two classes of theories have been examined 
depending on the relation between (w, fi, U) and (0, tp). Each of them has been 
studied without or with a perfect fiuid. 

Case A: without a perfect fiuid 

When the perfect fiuid is not present, we have for class 1 isotropisation: 
Case lA: i^{4>), A'(V') and U{(l),%li) 

A necessary condition for isotropisation of Bianchi type I model when two 
minimally and massive scalar fields are present will he that the two quanti- 
ties i^, = (j)U^U-'^{3 + 2w)-i/2 ^ ■4>U^U~'^{3 + 2^)-i/2 ^g„(^ ^0 ^^^g 
constants such as < 3- When isotropisation occurs and one of the two 

constants is non vanishing, the power law <i>i'^ i-i' is a late time attractor 
of the metric functions and the potential vanishes as t^^ . If the two constants 
vanish, the de Sitter Universe represents the late time attractor and the potential 
tends to a constant. 

If we put = strictly, we recover the results got in presence of a single 
scalar field without a perfect fiuid pi]. Results of case lA can be generalised 
for n scalar fields 4>i when their associated Brans-Dicke coupling functions u)i 
respectively depend on only 4>i (see annexe 1). For that, it is sufiicient to replace 
~^ Si^?- In the literature, it has been shown that the presence of 

multiple scalar fields could help to generate infiation: it is the assisted infla- 
tionjJlj. Sometimes it is the opposite which happens: the more scalar fields 
there are, the less likely the inflation occv,rs]2,^. It seems that it is this last be- 
haviour which arises for case lA: the more scalar fields there are, the more they 
contribute to the denominator of the power law to which the metric functions 
converge, and the less likely it will be larger than 1 and produce an accelerated 
expansion at late times. 

To evaluate the above conditions for a specific theory, we need to know the 
asymptotic behaviours for (j) and i/' such as we could calculate and ^^2- It 
comes: 
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asymptotic behaviours of the scalar fields for case lA 

The asymptotic behaviours of the two scalar fields when an isotropic state is 
reached are those of the two functions 4> and ip when fl — oo defined by: 

^ ^^^^^ (45) 



^ (3 + 2u)U 

Now we summarize our results concerning the second type of coupling for u), fi 
and U : 

Case 2A: U!{4>,i>), li{i>) and U{tp) 

There exists two equilibrium points Ei and E2 which may correspond to an 
isotropic equilibrium state with two minimally and massive scalar fields for the 
Bianchi type I model. The necessary conditions to reach equilibrium are ex- 
pressed with the two quantities = tpU^U~^{3 + 2^)~^/'^ and = 'il)U!^{3 + 
2w)-i(3 + 2^)-i/2; 

• For point El, it is necessary that < 3 and (3 — £^^51 — 2 ^ i^^l^^^dO. — > 
—00. When isotropisation occurs and if l^i,^ tends to a non vanishing con- 

stant, t "^i is the late time attractor of the metric functions and the poten- 
tial vanishes as t~^. If i^^ tends to zero, a de Sitter Universe is the late 
time attractor and the potential tends to a constant. If moreover ^i/,^ di- 
verges, an additional condition for isotropisation is £.,11.26^^^^ ^v-i)^ ^J^^^^^^ 
0. 

• For point E2, it is necessary that < 2£i/,2(^i/'i + 2i'^2)^^ < ^Vn + 
2£^2 7^ and £^^{£^.^ + 2i!^2) > 3. When isotropisation occurs and if 
£-Ki)i{£^i+2l^^)~^ tends to a non vanishing constant, the late time attractor 
of the metric functions is a power law of the proper time i(^*i+2'^»'2)(3^'*i) 
and the potential vanishes as t~^ . +2^^2)~^ vanishes, a de Sitter 
Universe is the late time attractor and the potential tends to a constant. 

Contrary to what happens with only a single scalar field, there are now two 
equilibirum points. The first one is such as the metric functions asymptotic 
behaviour only depends on ijj whereas for the second one, it depends on both 
scalar fields <f) and ip, the power law exponent being totally different from the 
previous case. The scalar fields asymptotic behaviours ahowing to calculate the 
quantities and are given in the following way: 

asymptotic behaviours of the scalar fields for case 2A 

The asymptotic behaviours of the two scalar fields when an isotropic state is 
reached are those of the two functions (j) and ip as Q. ^ —00 defined by: 

• for the equilibrium point Ei : 

^ = 12(j>{3 + 2a;)-i/2e(3-^'^i)"-'/^^i^'^^''^ (47) 
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• for the equilibrium point E2 : 




(49) 



[U{3 + 2uj)] 



U{3 + 2lj) 



e 



,6(0-J^o) 



(50) 



Let us examine the connection between our results and Wald's No Hair 
theorem Jl^. The latter states that initially expanding homogeneous models 
with a positive cosmological constant (except Bianchi type IX) and a stress en- 
ergy tensor satisfying the dominant and strong energy conditions, exponentially 
evolve to an isotropic de Sitter solution. The behaviour of Bianchi type IX 
model is similar if the cosmological constant is sufficiently large compared with 
spatial-curvature terms. Assuming that the Universe tends sufficiently fastly to 
its isotropic equilibrium state^, Wald's No Hair theorem can be generaHsed for 
the case 1 A to any form of potential and Brans-Dicke coupHng functions such as 

and tend to zero. For the case 2A and the equihbrium point Ei, only the 
vanishing of is necessary and, for the equilibrium point E2, the vanishing of 

(£^1 + 2i^^)~^ . For both cases, the potential tends to a constant, showing the 
stability of Wald theorem with respect to the presence of several scalar fields. 
Note that the relations between Bianchi models and Wald's No Hair theorem 
has been explored in the context of chaotic infiation in |27| . 

Case B: with a perfect fiuid 

When we take into account a perfect fluid, we get different conditions and met- 
ric functions asymptotic behaviours resulting from class 1 isotropisation. There 
exists two possible equilibrium points respectively corresponding to a vanishing 
or non vanishing k or equivalently the perfect fluid density parameter Q,m- For 
the flrst case, we have: 

Case IB: ui[4)), /i(V') (ind U{(f>,'ip)- 

A necessary condition for isotropisation of Bianchi type I model when 2 massive 
and minimally coupled scalar fields with a perfect fluid are considered and such 
as ri„i "> const ^ will he that the quantities = 4>U^U~^{?> + 2a;)~^/^ and 
— ■ipU^U^^{3 + 2fi)^^/'^ tend to some constants with £'^^ +^|^ > 3/2j. Then, 

when isotropisation occurs, the metric functions will always tend to t^-< and the 
potential will vanish as t^^ . When isotropisation arises such as Vim ~* 0, we 
recover the results of case lA (including the scalar fields asymptotic behaviours) 
hut the condition on -I- i?^^ is cast into £'^_^ + < 8/27. 

Hence, when fl„i const ^ 0, the metric functions asymptotic behaviour 
is the same as in presence of a single scalar fleld^^l- We flnd for the scalar 
flelds asymptotic behaviours: 

Asymptotic behaviours of the scalar flelds for case IB {k ^ 0) 

The asymptotic behaviours of the two scalar fields when an isotropic state is 

^i.e. considering that the assumptions on /(</>, 1/1) and {y, z, w, k) we explained in the section 
|3]are checked 
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reached with Ctm — * const ^ are those of the two functions (j) and i/} as 
f2 — > — oo defined by: 

: (3 + 2ti)c^^UU^ 

^ ^^(3 + 2m)</.2C/| + (3 + 2u;)V2/72 ^ 

; ^ o (3 + 2lu)<p^UU^ 

If we consider the second type of coupling, we get the following results: 
Case 2B: a;(0,V'), A*(V') and U{ij)). 

Let he the quantities l-^,, = V^f/^C/^HS + and = ii'cc;,.(3 + 2a;)-i(3 + 

2;u)~^/^. Some necessary conditions for isotropisation of Bianchi type I model 
when 2 massive and minimally coupled scalar fields with a perfect fluid are con- 
sidered and such as Q.m const ^ will he that tends to a constant with 

> 3/27 and (1 — 7/2)17 — 7 jC.^,2(.^^d^ —00 asVl^ —00. When isotropi- 

2 

sation arises, the metric functions will always tend to t '^- and the potential will 
vanish as t^"^ . When isotropisation arises such as flm 0, we recover the 
results of case 2 A (including the scalar fields asymptotic behaviours) but nec- 
essary reality conditions for isotropisation to Ei and E2 equilibrium points are 
cast into respectively < 3/27 and 1 — 7/2 < 2l^^{t^^^2i^^)~^ < 1- 

Once again when 0„i const 7^ 0, we recover the same behaviour for the 
metric functions as in the presence of a single scalar field despite the unusual 
form of w. To check the above limits and inequalities, again we need to know 
the scalar fields asymptotic behaviours. We have got: 

Asymptotic behaviours of the scalar fields for case 2B {k ^ 0) 
The asymptotic behaviours of the two scalar fields when an isotropic state is 
reached with Q.m const ^ are those of the two functions (j) and as 
f2 — > — 00 defined by: 

^ = (53) 

. ^ 120 ^3[(i-7/2)-7/W^i^rfn] (54^ 
^ v'3 + 2a; ^ ^ 



6 Applications 

To illustrate our calculation, we look for isotropisation conditions of some im- 
portant theories extensively studied in the literature. Remember that we have 
assumed a positive potential, the respect of the weak energy condition and 

7 G [1,2]. Then, in the following applications, when we will write that isotropi- 
sation is impossible, we must keep in mind that it could be wrong if one of 
the above assumptions were violated. These appHcations will be illustrated 
with numerical simulations done with a Runge-Kutta algorithm (order 5) im- 
plemented in Java. Java program and its user manual may be downloaded at 
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http:/ /luth2. obspm.fr/~etu/fay /stephaiie.htmH allows to integrate any hy- 
perextended scalar tensor theories (with varying G, w and U) with a perfect 
fluid for any class A Bianchi models written with the Lagrangian or Hamil- 
tonian(with the variables of this work) field equations. The equations system 
(I24I27|I is also implemented and any new equations system or numerical meth- 
ods may be easily added and should work with all the codes already written if 
user manual recommendations are followed. 

6.1 Hybrid inflation 

In the introduction we have connected the presence of two scalar fields with 
higher order theories or hybrid infiation. Hybrid infiation is studied in ^ with 
a scalar tensor theory defined by: 

(3 + 2w)0-2 ^ 2 (55) 
(3 + 2^)V"^ = 2 (56) 
U = 1 /4A(V'^ - A/2) + 1 /2m^(t? + 1 /2A (57) 



m, M , A and A' being some constants. It thus corresponds to cases lA and 
IB defined above. The same type of theory is also used in |^ for similar 
reasons and from the point of view of topological defects. The potential lf57jl 
has the symmetry (j) ^ —0 and 7/1 ^ —-0 and is the most general form of a 
renormalisable potential with this property, apart from the absence of a A"(/)^ 
term. For a flat FLRW model, inflation stops when the true vacuum state, 
which corresponds to a global minimum for the potential with {(f>,ip) — (0, M), 
is reached. When no perfect fluid is present, we calculate that and i^i,^ are 
respectively proportional to (j) and ip and write: 



2V20(to2 + A'V^^ 



A(M2 - 7/;2)2 + 202 (m2 + A'V'2 



(58) 



(59) 



A(M2 -i/;2)2 + 202(^7^2 + ;^/^2) 

Obviously, with {(f), ^p) = (0, M), we have — > and — -02 — > 0. Then, if we 
assume that the vanishing of is slower, faster or of the same order as M'^ — -02, 
we respectively flnd that i^-^, or the couple diverge. Then it is 

the same for the derivatives of the scalar flelds. The flrst graph of flgure ITTl 
represents a numerical integration of the scalar fields and illustrates this fact. 
Consequently, the couple (0, -0) = (0, M) represents an asymptotic state of true 
vacuum which can not occur with isotropisation of the Bianchi type / model. 
Moreover, numerical simulations show that the scalar fields are not defined as 
Q, ^ —00, thus confirming that this theory can not lead to isotropisation of 
the Universe. Such a result is interesting because early time infiation is often 
used to solve some problems of the standard big bang model such as the fiatness 
problem or the isotropy of the cosmological microwave background. However we 
see that starting from an anisotropic model, the hybrid infiation for the theory 
defined by Ij55l57|l is not able to isotropise the Universe. 

When a perfect fiuid is present, numerical simulations (second and third graphs 
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of figure indicate that 4> would oscillate to whereas %p would tend to a 
constant Mq different from M as SI ^ —oo. Thus, the potential tends to a 
constant and not to V~'' . Consequently isotropisation does not occur when 
k ^Q. Since it can not arise either when no perfect fluid is present, we conclude 
that, even when fc — » 0, isotropisation does not take place at late time. 
Hence class 1 isotropisation seems impossible for the theory of this section. 
Numerical simulations for the system Ij24l27ll confirm the result and do not 
show class 2 or 3 isotropisation either. 

6.2 High-order theories and compact ificat ion 

Another theory can be defined by the same forms of Brans-Dicke coupHng func- 
tions but with another form of potential: 

[/ — [/j^g--\/273fe0g-5y3/6fcVj^gy3/2-!/J _ ^gg^ 

with fc > and to > 0^. Such potentials appear when we compactify the 
space-time and cast high-order theories of gravity into relativistic forms. Hence 
in [3, conformal transformations are applied to the theory defined by 5* = 

/ (f' x\/Gz{^^R5^- aM^^ Ri^) and lead to the scalar tensor theory defined above 

with TO — 4/3, whereas if we consider the action S = J d^x-jGl{^^R. 



1.5 -t- 

we get a scalar tensor theory with to = 2. These actions 
are related to M-theory compactification. When no perfect fiuid is present, 
using and lUfill . we find that near isotropy: 



(61) 



-{2V3TOln fe^'''/2(5fc - 3to) - 5fc 



5(5fc — 3to) 

-(5fc-3TO)V'} (62) 



Since we consider /c > 0, -0 does not diverge to —00 otherwise the left mem- 
ber of equation ll62|l would be complex. Numerical simulations show that V 
tends to -l-oo when Q. —00 and then, we deduce from Ij62|l that -i/i 
-{5k - 3m){2V3)~^n. This limit will arise in ^ —00 if in the same time 
5k — 3m > 0. An illustration of the two scalar fields asymptotical behaviours has 
been plotted on the fourth graph of figure ITTl We calculate that the quantities 
and respectively tend to the constants —k/^/S and (3to — 5fc)/(2\/6). 
The necessary condition for isotropisation is thus (llfc^ — IOAjto + 3m?)/8 < 3. 
Assuming that (k,m) 7^ (0,0), the late time attractor of the metric functions 
is a power law of the proper time t^'*[**^ +(5A:-3m) ] ^ Hence, after some con- 
formal transformations, these theories derived from particle physics can lead to 
isotropisation of Bianchi type / model as illustrated by figure El 
When a perfect fiuid is present, numerical analysis of lf52ll shows that scalar 
fields are defined when 17 — > — 00 and may diverge. From the forms of cf) and 
V', it is easy to see that can not tend to —00 for positive k when Vl —00. 

^ These assumptions allow to simplify the study. 
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Figure 6: These figures, with — H in abscissa, represent successively the behaviours of 
{x,y,z,w,<p,ip,e^-^) for initial condition (x,y, z,w, (l>,il>) = (-0.49,0.25,-0.12,-0.15,0.14,0.23) 
and parameters (C/o,A:, m) — (3.2,1.25,-0.36) and a dust fluid. Note that is a constant 
— fc/\/3 — —0.721688. The last figure shows the vanishing of a, /3 and 7 derivatives with respect to 
the proper time as it should be in case of metric functions convergence to a power law of time. If 
we take m — —2.36, (llk^ — 10km + 3m^)/8 > 3 and class 1 isotropisation does not occurs since x 
tends to a non vanishing constant 



When Tp +00, it comes + — > (llfc^ — 10km + 3m^)/8 and thus this 
theory may isotropise to an equihbrium state whose nature depends on the value 
of this constant with respect to 3/2j. This case is illustrated on figuredwhere a 
numerical integration has been performed with (llfc^ — lOfcm + 3m^)/8 > 8/27. 
Numerical integration for the scalar fields theoretical asymptotical behaviours 
has also been plotted on the fifth graph of figure It also produces some so- 
lutions for which ip vanishes and (j) tends to a non vanishing constant but then, 
^^1 + ^'tpi diverges and class 1 isotropisation should not occur. 

6.3 A common quadratic potential for a complex scalar 
field 

The theories corresponding to cases 2A and 2B are related to the presence of 
complex scalar fields whose Lagrangian is most of times written as [201 1^ \29\ : 

L = i? + g'"'C,;C,.-^(ICn+im (63) 

By redefining the scalar field C as C = ip{V2m)e~'^"^'^ , it becomes: 

L^R+ l/25^''(V'2,/),^0,, + m- V.mV-,.) - C/(^') + (64) 

which corresponds to 3/2 + fJ- ~ l/2m~^V^ and 3/2 + co = 1/20^-0^. Since 
the potential depends on ip'^, its most simple and maybe natural form seems 
to be [/ = CC* = ip"^- It is often used in the literature for instance for scalar 
fields quantization in [201 or to study the genericity of infiation for spatially 
closed FLRW models in [22|. If we suppose that there is no perfect fiuid, then, 
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Fi gure 7: These figures, with — n in abscissa, represent successively the behaviours of 
{x, y,z,w, <p, ^, k, £^ J for initial condition (x,y, z,w, 0, -0) = (-0.49, 0.25, -0.12, -0.15, 0.14, 0.23) 
and parameters (t/o,/e,m) — (3.2,1.25,-2.36) and a dust fluid. Note that is a constant 
— fc/\/3 — —0.721688. As previously, last figure shows the vanishing of a, j3 and 7 derivatives 
with respect to the proper time. 



for El equilibrium point, we get V ±.2myj2{n — -^o): it is complex when 
— > —00 whereas, by definition, it should be real. For E2 equilibrium point, 
we get 7/; Voe"^^^^ whereas now (f> tends to a complex value instead of a 
real one. Consequently, for the theory defined by II64II with U — ijj^ , class 1 
isotropisation does not occur at late times. However, numerical simulations of 
equations H24I25II reveal that Universe may ("may" and not "must" since x ^ Q 
is a necessary but not sufficient condition for isotropisation.) undergoes a class 
3 isotropisation as shows on figure |HI with the characteristics explained at the 
beginning of this work. 

If now we assume that a perfect fluid is present, ip ^1"^!^ and l^^ diverges as 
g-3/27a. ^jjgjj class 1 isotropisation is not possible if fc 7^ 0. However, once again 
a class 3 isotropisation is possible with k oscillating to a constant as plotted on 
figure El If fc ^ 0, as shown above, class 1 isotropisation is impossible but not 
class 3. 

6.4 Topological defects 

Another type of potential has been used in |0] to study the formation of topo- 
logical defects after early time inflation. Its form is J7 = A/2(7/)^ — ?7^)^, i.e. 
the so-called wine bottle potential, with A and 77 some constants. If we as- 
sume that there is no perfect fluid, we calculate for E\ equihbrium point that 
i/)^ -7]'^ ProductLog{-'i]^'^e^^^"^ ^ (^^-^o))^ 0^ being an integration con- 
stant^. But this last quantity is negative when fl —00 and then, again, 
ip is asymptotically complex, which does not flt with its deflnition as a real 

^ ProducLog{z) gives the principal solution for w in z = we™ . 
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Figure 8: These figures, with —£7 in abscissa, represent successively the be- 
haviours of {x, x,y, z,w, (p,^, i^^, i^^) for initial condition {x^y,z,w,(p,ip) — 
(-0.70,0.25,-0.12,-0.15,0.14,0.23) and parameters m = -2.3. x is the only variable to 
reach equilibrium whereas y, z and w oscillates more and more as —Q increases. The scalar fields 
undergo damped oscillations whereas the oscillations for and increase. The last figure 
shows the vanishing of a, /3 and 7 derivatives with respect to the proper time as needed for 
isotropisation. Note that they oscillate. 

scalar field. For point E2, we also find that ip is complex as ^ —00 but if the 
integration constant is complex too. So for both Ei and E2 points, an isotropic 
equilibrium state of class 1 type can not be reached because at least one of the 
scalar fields is complex at late time. 

If now we assume that a perfect fiuid is present, we have — > ^'^/^li^-^^o} j^rf'^ 
with VIq an integration constant. Hence, diverges and class 1 isotropisation 
does not occur for the same reasons as in the previous application. 
However, once again, we have observed class 3 isotropisation with and without 
matter. In the first case, k tends to a constant with damped oscillations and we 
have observed that x but also z and the scalar fields could reach equilibrium. 
This is depicted on figure The same remarks apply to the second case. 
Overall the behaviours of the functions are the same as these shown on figures 

IHI 
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Fi gUrG 9l If we take into account a perfect fluid, k may reach a constant value during isotropisa- 
tion. 

6.5 Bose-Einstein condensate 

In [22], Bose-Einstein condensate is studied" with a potential of the form a^"^ + 
f3tjj'^. Again, assuming no perfect fluid, V is complex for Ei equilibrium point. 
Indeed, ip [a{2l3^^)] ^'"^ {ProductLog{a-^e^+^^"'^^°'''^'^~'f"''>)-l)^/^ with V'o 
an integration constant. Thus, when Q —00, the second square root is real 
only if a/3~^ < but then the first one is complex. For E2 equilibrium point, 
■0^ tends to the constant —al3~^ with a < and /3 > 0. In the same time, 
— 2(— 3/3Q!~^)^/^r2 + 4>o, 4>o being an integration constant. Calculating (^-^ 
and we get respectively that diverges and i^^ ~^ ±m\/— /3a~^. Hence, 
21^^ {£^^ + 2£^^)'^^ and y ^ Q. We could have a class 2 isotropisation 
although numerical simulations have failed to confirm it. 

If now we consider a perfect fiuid, we find ijP' — a(2/3)^^ ± (2/3)^^ (a^ + 
4/3e~'^'''*^^°~^^)^/^, r^o being an integration constant. Then, diverges and 
an isotropic stable state may be reached only if A; 0. From what we have 
found above, isotropisation could only occur for E2 point. However, since the 
vanishing of k needs 1 — 7/2 < 2£^2 (^'^^+2^^^ and the right member of this 
inequality is vanishing, we conclude that the theory should not undergo class 1 
isotropisation. 

Once again numerical simulations show class 3 isotropisation, with or without 
a perfect fiuid, and with the same behaviours as those shown on figures |H1 

We observe that all the theories dealing with complex scalar fields seem to 
reach isotropisation via class 3 mainly, whereas the others may reach it via class 
1. 

7 Conclusion 

We have studied the isotropisation of the fiat homogeneous Bianchi type / model 
filled with a perfect fiuid and two real scalar fields. This is an important issue 
because, as explained in sectional such theories are used to describe hybrid in- 
fiation, compactification mechanisms, topological defects or Bose-Einstein con- 
densate which may be related to primordial Universe. Taking the point of view 
that early Universe is anisotropic, the Lagrangian describing these theories have 
to be constrained to explain why isotropy arises and what looks like the Uni- 
verse isotropic state. 

To reach this goal, we have made the following assumptions: 

• We consider the scalar fields either such as and U{(l),ip) or 

■'The Lagrangian is different from JHSJ- 
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Fi gUrG lOl These figures, with — in abscissa, represent successively the be- 
haviours of (a:, ±, y, i, uj, 0, i/j, ' ) initial condition (x, y , 5:, u;, 0, = 
(0.49,0.25,-0.12,-0.15,0.14,0.23) and parameters (A, rj) = (0.25,0.25). x, z and the 
scalar fields reach equilibrium whereas undergoes non damped oscillations. The last figure 
shows the vanishing of a, fS and 7 derivatives with respect to the proper time. 



uj{(j),ip), iJ.{tp) and U{ip) since we thus recover a large number of theo- 
ries with two real scalar fields or one complex scalar field studied in the 
literature. 

• The weak energy condition is satisfied. 

• The potential, which may be considered as a variable cosmological con- 
stant, is positive. 

• Asymptoticalljf the density parameter of the scalar field should tend to 
a non vanishing constant value and the ratio of its pressure and energy 
density to a negative value in accordance with WMAP data. 

• The isotropic state is reached sufficiently fast. 

We have then found some necessary conditions for the Universe isotropisation 

under the form of some hmits and inequalies expressing with respect to the 
functions , i^^ ) ^i^i ^-^d of the scalar fields (j) and ip- The natural outcome 
of the Universe isotropisation may be described as 

• A De sitter Universe with a non vanishing cosmological constant 
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• An Einstein - De Sitter Universe (e ^ ^ t^i) with a vanishing cosmologi- 
cal constant (J7 t^^) and a non vanishing perfect fluid density parameter 

- 

• A power law expanding Universe (e~^ t™, with m the limit of a deter- 
mined function of the scalar fields) with a vanishing cosmological constant 

and a vanishing perfect fiuid density parameter flm- 

Note that the potential always tends to a constant or decreases as whatever 
the forms of uj, fi and U. In this last case, if the Universe is 15 Gys old, the 
cosmological constant should be 4.96.10~^^cto~^ in agreement with super- 
novae observations. 

When there is no perfect fluid and (^(0), and U{4>,ip), the results 

generalise those of |1SJ, where a single scalar field is considered, for any num- 
ber of minimally coupled scalar fields whose associated Brans-Dicke coupling 
functions uji only depend on 0^ : a necessary condition for isotropisation is that 
Yli^^i tends to a constant smaller than 3. If the constant is vanishing, the 
Universe tends to a De Sitter model otherwise the metric functions increase 
as t ' 'I'i . When uj{(f>,ip), l^{ip) and U{ip), the results are different because 
now the factor in front of the 4> field kinetic term contains the ip field. Hence, 
we find two equilibrium points and the power laws representing the asymptotic 
behaviours of the metric functions when isotropisation occurs are different from 
the previous case or what we had found in {15j. 

Considering a perfect fluid modifies the necessary conditions for isotropy 
even when its density parameter Sim tends to vanish. However, in this last case, 
the asymptotical behaviours of the metric functions and potential are the same 
as without it whereas if Qm tends to a nonvanishing constant, the metric func- 

2 

tions behave as if they were no scalar field at all, i.e as t^^ . Hence, when fi^ 
and rim are asymptotically of the same order, the expansion is decelerated thus 
preventing to solve the coincidence problem^. 

From an observational point of view, this paper shows that the presence of 
several minimally coupled scalar fields would not be detectable by dynamical 
observations of a nearly isotropic Universe since the dynamical behaviours of 
the metric functions or potential are of the same nature as in the presence of a 
single one, thus showing a degeneracy problem. 

We have applied our results to several theories and shown, considering the 
above assumptions, that the model of hybrid infiation considered in [2] does 
not lead the Universe to an isotropic state on the contrary to some theories 
coming from compactification process and studied in [Jj. All the theories with 
a complex scalar field and related to scalar fields quantization|2D], topological 
defects(2Ii| or Bose-Eisntein condensate[22 do not undergo a class 1 isotropisa- 
tion but a class 3, showing among others strongly oscillating behaviours of one 
or two scalar fields or even of the perfect fiuid density parameter. 
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Figure 11: Some scalar fields numerical integrations for the following applica- 
tions of section 4: hybrid infiation (m = 1, M = 1, A = 1, A' = 5, = 1, 
■^(l) = 1) without and with a perfect fluid, high-order theories {m — I, k = 2, 
0(1) = —2.5, ip{l) = 0.70) without and with a perfect fluid. 
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A Generalisation of case lA for n scalar fields 

If we consider the presence of n scalar flelds 0„ and we use the following vari- 
ables: 

X = H-^ (65) 



y 



= Ve-«2[/i/-i (66) 
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Zi = P4>M^ + 2iJi)-^^^H-^ (67) 

i varying from 1 to n, we get the following first order equations system from the 
Hamiltonian equations: 

X = 2,R^y^x (68) 

n 

y = 3y{2j2£iZi + RV-l) (69) 



Zi = y^R\3zi - l/2li) + 12 5^ 0^ " 12 ^H^j^i (70) 

with li = (t)iU4.,U-\3 + 2wj)-i/2 and £ij = (t>iOJj4>,{2> + 2uJi)-^/^{2, + 2ujj)-^. 
If we assume that each Brans-Dicke coupHng function Wi only depends on the 
scalar field as for the case 1, Cij = when i ^ j and the equilibrium points 
are {x,y,zi, ..,Zn) = (0, ±(3 - EiU )'^'(v^'R)"'' V6^i, -, 1/6C). K is thus 
possible to generalise the results of case 1 by replacing + by Xl"=i ^1- 



B With a perfect fluid 

Equilibrium points calculus when l^^ = 







The equihbrium points are defined by the following (a;,y, values: 

• ^1 = (0,0,0,0) 

• E2,3 = (0,±[9(2 - 7)7 



^ti + 12(7 



2 (37 



2^il)]'/'[4l+4l]"'^'(2^/3i^)-^f0l(6-37+2^L+24l) fl2(^L + " ) 



£^i(6-37 + 2£2^ + 2£2^) 12(£2^+€2J ) 



1/2 



. ^;4,5 = (0,±l/2^/3i^-l [7(2-7)(^|,+^2^)-iJ , 1/47^(^2^ 
El may correspond to a class 2 isotropisation. For the other equihbrium points, 



the constraint implies k 



1 



, 37 



that is real and non vanishing if 



> 3/2^- ^2,3 points will be real only if G [3/2(7 - 2), 3/27] but 

this condition is incompatible with above constraint on k. Consequently they 
are eliminated from further considerations. £^4,5 points are real since 7 e [1,2[ 
and thus will be the only one we will consider. 

Equihbrium points calculus when = i^^ = 

We find 11 equilibrium points that we introduce in the constraint equation 
to determine the form for k. They can be divided into 3 groups: 



• First group: 

The constraint requires k^ 



1 and then equilibrium points are given by: 
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- Ei = (0,0,0,0) 

- S6.7 = (0,0,±l/ 



-(7-2) 



2^-2 
i/)2 



1/2 



,1/8(2-7)^-2) 



The El point is similar to that of the previous section and we make the 
same remarks. EquiHbrium points Eqj are complex and thus eliminated 
from further considerations. 

• Second group: 

The constraint requires 1 ~ 3/2^£^'^. Introducing this value in the 

equilibrium points, we get: 



4/4 



1/2 



,0, 



- = (0, ±l/2i?-i/^V37(2 - 7), 0, 1/47^^1'). 

- Ei,5 = (0, ±(2V3i?/^i)-i [97(2 - 7) + 12(7 - 1)£2 ^ 
(12/^i)-H6-37 + 2/2^)) 

i?2,3 points are real for the considered values of 7. i?4,5 points are real 



if 3/27 e 



which is not compatible with a real k arising for 



3/27 < /^j. Hence, £"4,5 points are eliminated. 
• Third group 

In this last group, the constraint requires fc = or fc ^ 0. Then equilib- 
rium points values and x asymptotic behaviour are the same as in section 
14.11 although isotropisation conditions are modified as shown in subsection 

UTS 

It follows that only i?2,3 equilibrium points may represent an isotropic stable 
state when i"^-^ > 3/27 and k is strictly or asymptotically different from zero. 
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